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Abstract. Refractory black carbon aerosols (rBC) from 
biomass burning and fossil fuel combustion are deposited 
to the Antarctic ice sheet and preserve a history of emis- 
sions and long-range transport from low- and mid-latitudes. 
Antarctic ice core rBC records may thus provide informa- 
tion with respect to past combustion aerosol emissions and 
atmospheric circulation. Here, we present six East Antarc- 
tic ice core records of rBC concentrations and fluxes cov- 
ering the last two centuries with approximately annual res- 
olution (cal. yr. 1800 to 2000). The ice cores were drilled 
in disparate regions of the high East Antarctic ice sheet, 
at different elevations and net snow accumulation rates. 
Annual rBC concentrations were log-normally distributed 
and geometric means of annual concentrations ranged from 
0.10 to 0.18pgkg _1 . Average rBC fluxes were deter- 
mined over the time periods 1800 to 2000 and 1963 to 
2000 and ranged from 3.4 to 15.5 pgm -2 a -1 and 3.6 to 
21.8 pgm -2 a -1 , respectively. Geometric mean concentra- 
tions spanning 1800 to 2000 increased linearly with eleva- 
tion at a rate of 0.025 pg kg - 1 /500 m. Spectral analysis of the 
records revealed significant decadal-scale variability, which 
at several sites was comparable to decadal ENSO variability. 


1 Introduction 

Nanoparticles of refractory black carbon (rBC, soot) aerosols 
are emitted to the atmosphere during fires and fossil fuel 
combustion and transported over long distances at the 
hemisphere-scale (Seiler and Crutzen, 1980; Andreae et ah, 
2005; Crutzen and Andreae, 1990). Because of their strong 
light absorption properties, rBC nanoparticles alter air tem- 
perature, snow albedo and impact climate (Moosmiiller et 
ah, 2009; Ramanathan and Carmichael, 2008; Ramanathan 
et ah, 2001; Jacobson, 2001; Flanner et ah, 2007; Penner et 
ah, 2002). The incorporation of future emissions scenarios is 
therefore a key parameter for global climate modellers (Ra- 
manathan and Carmichael, 2008). 

However, even if surveys of rBC atmospheric concentra- 
tions and fire occurrence are being constructed for the last 
few decades thanks to the development of new satellite tools 
(Ito and Penner, 2004; Chung et ah, 2005), the history of rBC 
fire emissions over last two centuries is incomplete (Mouillot 
and Field, 2005). Indeed, a detailed fire history over periods 
covering the pre-industrial time to the modern era is lacking, 
notably because common fire proxies such as charcoal de- 
posits or gas emissions reconstructions from ice cores cannot 
reconstruct large spatial variability at inter-annual time scales 
(Marlon et ah, 2008; Wang et ah, 2010). Nonetheless, fire 
regimes are likely to have been modified by anthropogenic 
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activities such as tropical deforestation, land-clearing prac- 
tices or through natural changes in precipitation and tem- 
perature patterns with the changing climate (Dube, 2009; 
Nitschke and Innes, 2008; Bowman et al., 2009). 

Recently, several studies have directly used rBC from ice 
cores to estimate emissions from past combustion, spanning 
the transition from the pre -industrial to the modern era. In the 
Northern Hemisphere, ice cores from Greenland and the Hi- 
malayas have shown the impact of coal combustion on rBC 
snow concentrations (McConnell et al., 2007; McConnell, 
2010; Kaspari et al., 2011). In the Southern Hemisphere 
(SH), two high resolution rBC records from Antarctica have 
recently been used to investigate the evolution of forest and 
grass fires (Bisiaux et al., 2012). The records, from the West 
Antarctic Ice Sheet (WAIS) and Law Dome, spanned the time 
period 1850-2001 and showed large-scale changes in rBC 
deposition linked to climatic oscillations such as El Nino 
Southern Oscillation (ENSO) and the anthropogenic modi- 
fication of fire regimes. 

However, the contribution of each continent to Antarctic 
rBC deposition remains unknown, as well as the influence 
of atmospheric transport to sites located in the Atlantic sec- 
tor of Antarctica. Ice core sites located at high elevation on 
the East Antarctic Plateau in East Antarctica may record dif- 
ferent variability than the sites previously studied on coastal 
East Antarctica and the West Antarctic Ice Sheet (Bisiaux et 
al., 2012). The present study uses rBC paleorecords from the 
Atlantic sector to investigate spatial and temporal variability 
of rBC deposition and the link with emissions sources. We 
use six ice core records from Dronning Maud Land (Fig. 1) 
and focus on the cal. yr. period 1800-2000, covering pre- 
industrial and modern eras. These sites have much lower 
annual snow accumulation rates (20-60 kg m~ 2 a -1 ) com- 
pared to the sites investigated by Bisiaux et al. (2012) (150- 
200 kg m -2 a -1 ) and have a lower temporal resolution (an- 
nual to multi-annual). We use accumulation calculations 
from Anschutz et al. (2011) to estimate rBC fluxes at these 
sites and the importance of precipitation on rBC concentra- 
tions. Periodic oscillations in the rBC records are investi- 
gated through spectral and multiple regression analysis and 
compared to sodium (Na) records (measured simultaneously) 
to evaluate the influence of transport in the observed rBC 
variability. 

2 Drilling sites and methods 

2.1 Drilling site locations and characteristics 

The sites are all located on the East Antarctic Plateau (el- 
evation > 2500 m), in Dronning Maud Land (Fig. 1). The 
cores were drilled during two exploratory traverses by a 
Norwegian-American team (NUS) in the summer of 2007 
(cores “07-X”) and in the summer of 2008 (cores “08-X’). 
Latitudes, longitudes and elevations are compiled in Table 1 . 


0°0'0" 1 5°0'0"E 



105°0’0"W 105°0'0"E 


Fig. 1. Map of the traverse route 2007/2008 (black line) and 
2008/2009 (blue line) with drill described in this study sites from 
both legs indicated (NUS07-X and NUS08-X). Relevant stations in 
the area of investigation are shown as well. Elevation contour lines 
are in 100 m intervals. The map was compiled by K. Langley and 
S. Tronstad (Norwegian Polar Institute) and adapted by Anschutz et 
al. (201 1) for this study. The locations of the WAIS and Law Dome 
drilling sites are shown on the inset as red dots. 


Cores NUS08-4 and NUS08-5 were drilled in adjacent loca- 
tions and only ~17 km apart. All cores are firn cores of total 
length under 90 m, of which we present the top part, down 
to a depth corresponding to cal. yr. 1800 (Table 1). How- 
ever, the top dates are not common to all cores (due largely 
to fragile near-surface firn sections) and range from cal. yr. 
1989 for NUS 07-5 to cal. yr. 2008 for NUS 07-7(Table 1). 

2.2 Ice core analysis 

Longitudinal sections of the NUS ice cores were analysed 
from 2008 to 2010 at the Desert Research Institute on an ice- 
core melter continuous flow analysis system with in-line rBC 
measurements. The inner 1 cm 2 of the sections were used for 
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Table 1 . Sites coordinates. l'BC concentrations and fluxes, annual accumulation for this study and previous study by Bisiaux et al. (2012), 
for two periods of time (since 1800 and since 1963). 


This study | Bisiaux et al. (2012) 




NUS07-1 

NUS07-2 

NUS07-5 

NUS07-7 

NUS08-4 

NUS08-5 

WAIS (WDC06A) 

Law Dome (DSSW19K) 

Lat - long 


73°43 / S- 

76°04 / S- 

78°39' S- 

82°04' S- 

82°49' S- 

82° 3 8' S- 

79°46' S- 

66°73' S- 



07°59' E 

22°28' E 

35°38 / E 

54°53'E 

18°54'E 

17°52'E 

112°08 / W 

112°83'E 

Elevation (in m a.s.l.) 


3174 

3582 

3619 

3725 

2552 

2544 

1766 

1390 

Depth of yr 1800 for this 


22.2/30.5 

16.4/90.4 

12.6/89.5 

14.6/90.6 

18.05/30 

17.5/92 

- 

- 

Study/total core depth (in m) 










Period covered 


1800-2006 

1800-1993 

1800-1989 

1800-2008 

1800-2004 

1800-1993 

1850-2001 

1850-2001 

Number of data points 


2154 

1501 

1111 

1308 

1672 

1616 

4860 

2883 

Annual rBC cone, (geometric 

...since 1800 

0.16 

0.12 

0.14 

0.18 

0.10 

0.11 

0.08 

0.09 

mean) in pg kg - 1 . . . 

. . . range (2<r) a 

0.09 to 0.26 

0.07 to 0.19 

0.08 to 0.26 

0.12 to 0.27 

0.06 to 0.18 

0.07 to 0.18 

0.05 to 0.12 

0.05 to 0.2 


. . . since 1963 

0.14 

0.14 

0.18 

0.19 

0.15 

0.12 

0.08 

0.07 


. . . range (2cr) a 

0.08 to 0.27 

0.08 to 0.24 

0.14 to 0.24 

0.13 to 0.29 

0.08 to 0.26 

0.08 to 0.20 

0.05 to 0.12 

0.04 to 0.15 

Annual accumulation* 5 

...since 1815 

52.0 ± 2.0 

33.0 ± 0.7 

24.0 ± 0.5 

29.4 ± 0.6 

36.7 ± 0.9 

35.0 ± 0.8 

200 ±3.4 

150 ±3.1 

in kgm -2 a - *. . . 

. . . since 1963 

55.9 ± 3.9 

28.0 ± 2.0 

20.1 ± 1.4 

26.1 ± 1.9 

36.1 ± 2.1 

37.6 ± 2.3 

- 

- 

Annual rBC fluxes (geometric 

...since 1800c 

8.3 

3.9 

3.4 

5.3 

3.7 

3.9 

16 

13.5 

mean) in pg m -2 a - * . . . 

. . . range (2 tr) 

4.6 to 14.2 

2.5 to 6.2 

1.8 to 6.3 

3.5 to 8.0 

2.1 to 6.9 

2.2 to 6.5 

9.8 to 24.4 

7.3 to 30.6 


...since 1963 

7.8 

3.9 

3.62 

5 

5.4 

4.5 

- 

- 


. . . range (2 tr) 

4.0 to 15.8 

2.4 to 7.3 

2.6 to 5.2 

3.1 to 8.0 

2.7 to 10.0 

2.7 to 8.0 

- 

- 


a Multiplicative standard deviation. 
b From Anschutz et al. 2011. 

c We assume same accumulation rate for the 1800-1815 periods. 


the rBC analysis according to the method previously used 
in McConnell et al. (2007) and McConnell (2010), and de- 
scribed in detail by Bisiaux et al. (2012). The rBC analy- 
sis consisted of an ultrasonic nebulizer/desolvation system 
(CETAC UT5000) coupled to a Single Particle Soot Pho- 
tometer (SP2, Droplet Measurement Technologies, Boulder, 
Colorado). In this system, ice core meltwater was nebulized 
and desolvated to form a dry aerosol. The aerosol then passed 
into the SP2, where individual particles in the diameter range 
''-70-400 mn were heated up to incandescence by an Nd- 
YAG laser (1064 nm), and the emitted radiation was mea- 
sured by optical detectors (photomultiplier tubes). Individual 
particle masses were determined using calibration data gen- 
erated from the introduction of rBC particles of known mass 
directly into the SP2. Additional calibrations of the SP2 cou- 
pled to the ultrasonic nebulizer (USN) were performed daily 
using rBC colloids. These calibrations were used to account 
for rBC losses in the USN. Depth resolution of the analytical 
system was estimated at 1 cm. 

2.3 Ice core dating 

Annual layer counting was not possible for these NUS cores 
and dating was based on the identification and mapping be- 
tween ice cores of a number of chemical markers correspond- 
ing to explosive volcanic eruptions. Specifically, we used 
continuous, high-depth-resolution measurements of non-sea- 
salt sulphur (nssS) in the WAIS Divide (WDC06A) and NUS 
ice cores to identify layers with significant volcanic sulphur 
concentrations. Dating of the volcanic horizons was per- 
formed using annual layer counting in the WAIS Divide ice 
core (Bisiaux et al., 2012). The WAIS core volcanic chronol- 
ogy was then applied to the NUS cores’ volcanic horizons 


(Table 1). Although well-known large volcanic events (e.g. 
Tambora, Krakatoa, Agung) were included in the chronol- 
ogy, smaller, less-known volcanic events were used to fill in 
the depth-age relationships where possible (a total of seven 
tie points: cal. years 1810, 1816, 1837, 1863/4, 1884, 1965, 
1991/2). Dating between the volcanic horizons assumed uni- 
form accumulation between horizons. Snow accumulation 
rates (net) at each site (Table 1) were based on these depth- 
age relationships and ice-core depth/density profiles (An- 
schutz et al., 2009, 2011) and restricted to average rates be- 
tween the dating horizons (which included the ice cap sur- 
face). 

The number of rBC data points varied from 6 to 10 per 
year. However, the dating uncertainty is likely on the or- 
der of several years due to the low annual snow accumula- 
tion rates and physical processes, such as the redistribution 
of snow by wind. Overall, the dating uncertainty may be as 
large as ±5 yr between the dating horizons. The depth-age 
relationship of the NUS 07-1 core benefited from previously 
published accumulation rate data (Isaksson et al., 1999) and 
was used as the basis of further refinement of the 08-5 and 
08-4 depth-age relationship using decadal trends in the rBC 
records. The refinement process consisted of warping the 
08-4 and 08-5 rBC non-linear trend (described in Sect. 2.4) 
peaks and troughs to the 07-1 non-linear rBC trend (within 
the constraints of the volcanic chronology) using the Anal- 
yseries software (Paillard et al., 1996) lineage tool. The re- 
fined depth-age relationships resulted in coherent non-linear 
trends between the 08-4 and 08-5 records, which were drilled 
in the same region (separated by 17 km). 


www.atmos-chem-phys.net/12/3799/2012/ 


Atmos. Chem. Phys., 12, 3799-3808, 2012 



3802 


M. M. Bisiaux et al.: Variability of black carbon deposition to the East Antarctic Plateau 


2.4 Data analysis 

2.4.1 Period studied 

The study focused on two time periods, including 1800 to 
1990-2003 (time period depending on the ice core consid- 
ered) and 1963 to 1990-2003. These time periods were se- 
lected to coincide with the snow accumulation calculations 
by Anschutz et al. (2011). For the purpose of this study, the 
snow accumulation rate from the time period 1809 to 1815 
was assumed to extend to 1800. 

2.4.2 Concentrations 

The frequency distributions of the ice core rBC concentra- 
tions were determined to be lognormal. Statistics, such as 
the average and standard deviation, and mathematical op- 
erations such as re-sampling, smoothing, trend and spec- 
tral analysis, etc. were therefore performed in log-space 
(Limpert et al., 2001) to give “robust” geometric rather than 
arithmetic statistics. Note that here the geometric standard 
deviation is the multiplicative standard deviation (Limpert 
et al., 2001), which covers 68.3% of the variability, cor- 
responding to ermin conc = geometric mean ■ geometric stan- 
dard deviation, and crmax conc = geometric mean/geometric 
standard deviation. Prior to spectral analysis the concentra- 
tion data were re-sampled to an annual time scale using the 
“robust” piecewise linear interpolation method described by 
Paillard et al. (1996). A 21-yr smoothing window was used 
to capture decadal-scale variability and was performed using 
the non-parametric Nadaraya- Watson kernel density regres- 
sion method (Watson, 1964; Nadaraya, 1965). The time win- 
dow used for the smoothing represents a trade-off between 
the temporal resolution of the different records and the time 
span of the records, and variability in the later half of the 20th 
century. Monotonic trends were estimated using the nonpara- 
metric Thiel-Sen approach (Gilbert, 1987; Onoz and Bayazit, 
2003), while non-linear trends were calculated with singu- 
lar spectrum analysis (SSA) using Kspectra software (Ghil et 
al., 2002) at a 95 % Kendall level of confidence. Significant 
non-linear trends were reconstructed using the Kspectra SSA 
tool, which partially reconstructs the time series based on a 
linear combination of trend principal components (Ghil et al., 
2002). The non-linear trend data were then back-transformed 
from log-space and normalized (Z-scores calculated using 
the arithmetic mean and standard deviation of the complete 
trend reconstruction) for easier comparison of variability be- 
tween records. 

2.4.3 Deposition fluxes 

Flux calculations were based on longer-term average snow 
accumulation estimates obtained by Anschutz et al. (2011) 
from the volcanic chronology. Atmospheric fluxes of 
rBC were estimated by multiplying annual rBC concen- 
trations and accumulation rates corresponding to the two 


time periods chosen. Fluxes uncertainties were esti- 
mated as erminflux = crmin C0n c • {accumulation — ermin accu ) 
and ermaxflux = ermax conc • (accumulation + crmax aC cu)- 

2.4.4 Spectral analysis 

Spectral analysis was conducted using Analyseries software 
(Paillard et al., 1996). Significant periodic oscillations in the 
rBC and Na ice core records and their spectral coherence 
were investigated using the Blackman-Tukey method with a 
Bartlett window. In this study we define coherence as the 
fraction of common variance between two time series x and 
y through a linear relation, considering non-zero when co- 
efficients reach values >0.38 (Paillard et al., 1996). Here, 
we use raw data re-sampled to a 0.4 yr step with piecewise 
linear interpolation to perform the calculations in order to re- 
move lower frequencies, but to keep the eventual annual cy- 
cles. Coherence coefficients are given with 3 levels of con- 
fidence (low, medium and high), of which we only present 
the medium level, and only for coefficients >0.38 in Fig. 6. 
Principal components #1 extracted from records are calcu- 
lated with an embedding dimension of 20 and theory from 
Vautard and Ghil (1989). 

3 Results and discussion 
3.1 Concentrations and fluxes 

Concentrations of rBC were found to be log-normally dis- 
tributed with geometric means of annual concentrations 
ranging from 0.10 to 0.18 pg kg -1 since 1800 and from 0.12 
to 0.19pgkg _1 since 1963(Table 1). Overall, the NUS ice 
core geometric mean rBC concentrations were higher than 
the concentrations previously determined at lower elevation 
sites for the same time periods at WAIS (WDC06A) and Law 
Dome (DSSW19K) (Bisiaux et al., 2012). We attribute this 
difference to the very low annual snow accumulation rates on 
the plateau (Table 1) and thus to a low rBC dilution by snow, 
suggesting a significant fraction of dry versus wet deposited 
rBC. However, annual rBC fluxes estimates since 1800 are 
still lower (3.4 to 8.3pgm _2 a _1 ) than those determined for 
the WDC06A and DSSW19K (Table 1). 

The time series of annual and 21-yr smoothed rBC con- 
centrations at all NUS sites are shown in Fig. 2. Signif- 
icant monotonic trends, with superimposed decadal vari- 
ability (for the whole period 1800 - onwards, Mann- 
Kendall test double-sided p-values <0.0001), were found 
in annual rBC concentrations at sites 07-5, 08-4 and 08-5 
(Fig. 3a). The trends at these sites represented an increase of 
~0.03 ±0.01pgkg _1 /100yr. Comparison of the records at 
annual resolution revealed no significant cross-correlations, 
including for sites 08-4 and 08-5, which were within 17 km 
of each other. We conclude that surface processes (accu- 
mulation rate variability, sastrugi, blowing snow, and snow 
sublimation) prevent the analysis of rBC signal at annual 
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Fig. 2. Time series of rBC concentrations. Black line is annual 
(piece-wise linear integration interpolation of log raw data) and red 
line is 21 yr k-smooth on annual (calculated in log space). The pe- 
riods of relatively low values (1890-1910) and high values (1920- 
1940), as described in Fig. 3, are indicated as shaded areas. 


time scales and restrain the resolution to interannual/decadal- 
scales (Pomeroy et al., 1999). 

Decadal-scale variability was investigated using singular 
spectrum analysis non-linear trend reconstruction (Ghil et 
al., 2002). Significant non-linear trends over the entire pe- 
riod (1800 - onwards, p < 0.05, Mann-Kendall trend test) 
are shown in Fig. 3b-c (normalized as Z-scores). Non- 
linear trends from sites 08-4 and 08-5 (which were indepen- 
dently mapped to the 07-1 time scale) were highly correlated 
(r = 0.64, r 2 = 0.41, n — 195, p < 0.01). Correlation coeffi- 
cients for the other records were insignificant, but with some 
common features. Comparison of these non-linear trends, cf. 
Fig. 3b-c, revealed a period of low concentrations from cal. 
yr. 1890 to ~1915, common to observations made at WAIS 
and Law Dome (Bisiaux et al., 2012). Here, however, this 
drop is followed by a period of relatively high concentra- 
tions until ~1940 and peaking locally in the 1930s. While 
this peak was detected in the high resolution record from 
Law Dome (Bisiaux et al., 2012), it was absent from the 
WAIS record. With the exception of site 07-2 (Fig. 3b-c), 



1790 1840 1890 1940 1990 

Cal. years 

Fig. 3. (a) Monotonic trends for sites 07-5, 08-4 and 08-5 (Kendall 
significance = 99 %). (b, c, d) Non-linear trends normalized as 
Z-scores (Kspectra software, Kendall significance = 95 %) for the 
six NUS rBC records, as a function of time. Corresponding frac- 
tion of record variability is indicated next to record name (%). (b) 
Comparison of twin sites 08-4 and 08-5, re-scaled from 07-1 dating 
(plain curve). Original dating is shown as dotted line, (c) Other 
three records 07-2, 07-5 and 07-7. Shaded areas highlight specif- 
ically common features and/or trends, (d) Comparison of trends 
(Z-scores) from sites 07-1, 08-4, 08-5 with NADA variance (ENSO 
long term variability) in dotted line, scale inverted. 


the NUS rBC records also lacked the period of low variance, 
from ~1940 to ~ 1980, found in the WAIS and Law Dome 
records. Finally, the last 20 yr (1980-2000) show an increas- 
ing trend for the cores’ recording this period, which was also 
noted by Bisiaux et al. (2012) for WAIS and Law Dome. 

3.2 Effect of elevation 

The geometric average of the annual rBC concentrations at 
each site was found to increase linearly with elevation, cf. Ta- 
ble 1 for site elevation. Linear correlation coefficient (r) was 
0.81 when concentrations were averaged since 1800 ( r 2 — 
0.67, n =8, p = 0.01), and 0.92 when concentrations were 
averaged since 1963 ( r 2 = 0.86, n — 8, p < 0.01) (Fig. 4a). 
The slope of this linear regression corresponds to an increase 
of 0.015 and 0.025 pg for an elevation gain of 500 m, respec- 
tively. 
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Fig. 4. Geometric mean rBC (a) and Na (b) concentrations ~since cal. yr. 1800 (top) and 1963 (bottom) as a function of site altitude. For Na, 
regression line was calculated without the maritime site “Law Dome”, (c) Geometric mean rBC concentrations as a function of accumulation 
after cal. yr. 1800 (top) and 1963 (bottom) (Anschutz et al., 2011). Uncertainty bars refer to Table 1. Regression lines and coefficients are 
based on geometric mean values, r 2 indicates the coefficient of determination. * Values for Law Dome and WAIS from Bisiaux et al. (2012). 


The increase of rBC with elevation was previously ob- 
served in the southern latitude atmospheres by Schwarz 
et al. (2010) and modelled for Arctic snow by Skeie et 
al. (201 1). Stohl (2006) also modelled increased atmospheric 
loading in the southern latitudes and attributed this increase 
to rBC transport from lower latitudes towards the ice cap, 
along isentropic trajectories that may not reach the surface of 
the region lower than the plateau and remain at higher eleva- 
tions. 

To investigate whether rBC transport to the plateau is mod- 
ulated by the intrusion of marine air masses, the variabil- 
ity of Na concentrations (geometric means) with elevation 
was also investigated. In this case, no significant relationship 
was found (r 2 — 0.05, n = 7 , p > 0.05)(Fig. 4b). This con- 
firms observations previously made by Bertler et al. (2005), 
showing no link between increased elevation and Na con- 
centrations for altitudes above 2000 m. We hypothesize that 
this absence of correlation with elevation for Na and pres- 
ence of correlation for rBC are due to both a difference in 
the sources of Na and rBC aerosols and to a difference in 
atmospheric transport. Indeed, the main sources of Na are 
marine aerosols, which are transported to the East Antarctic 
Plateau by low pressure systems (Sneed et al., 2011) and dry 
deposited (Fischer et al., 2007). Transport processes associ- 
ated with rBC, therefore, appear to be different from those 
associated with Na. This suggests that rBC inputs to the 
atmosphere of the East Antarctic Plateau are not controlled 
by the intrusion of marine air masses and that transport in 
the upper troposphere may be important. Vertical profiles of 


rBC in the near Antarctic atmosphere reported by Schwarz 
et al. (2010), who found that rBC increased with altitude. 
Here, wet removal processes limit the lifetime of rBC near 
the boundary layer, while dry air in the upper atmosphere 
increases the rBC residence time. 

Snow water accumulation rates, on the contrary, do show 
a significant inverse trend with rBC concentrations, but 
only from 1963 onwards (r 2 = 0.74, r — 0.85, n — 8, p < 
0.01)(Fig. 4c). However, this correlation is determined 
mainly by the WAIS and Law Dome data points, with NUS 
sites clustering around the same values (Fig. 4c). Uncertain- 
ties inherent in the net snow accumulation rate must also be 
considered. Acknowledging these caveats, the slope of the 
linear regression of 0.030 pg in rBC for a 50mm decrease in 
accumulation can be compared with the increase of 0.025 pg 
rBC estimated for every 500 m in elevation, for the time pe- 
riod from 1800 to the present (Fig. 4a, top). Thus, for the 
two time periods shown in Fig. 4, the change in elevation 
may explain ~80 % of the difference in rBC geometric mean 
concentrations. Therefore, we suggest that the main process 
controlling the spatial differences in geometric rBC snow 
concentrations between the sites of the Antarctic Plateau is 
the decrease in accumulation and corresponding increase in 
rBC dry deposition, inducing less dilution of the particles. 
This relationship may explain the monotonic trend found for 
the record 07-5 (Sect. 3.1), which exhibits a decrease in ac- 
cumulation rate from the period 1 8 1 5 - onwards to the period 
1963 - onwards (Anschutz et al., 2011; Isaksson et al., 1999), 
cf. Table 1 . However, it is not as clear for the two other sites 
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Fig. 5. Cross correlation coefficients between Na and rBC from the 
same record. Na leads rBC for delays >0 and rBC leads Na for 
delays <0. 


displaying significant increasing monotonic trends, 08-4 and 
08-5, but no strong trend in accumulations(Table 1). 

3.3 Influence of transport 

For the sites 07-1, 08-4 and 08-5, cross-correlations of the 
annual rBC and Na suggest common high frequency vari- 
ability between the rBC and Na species at each site without 
leads or lags (Fig. 5). These data suggest a transport compo- 
nent linked to some of the high frequency variability. Non- 
linear low frequency trends similar to those found in the rBC 
records (Fig. 3b-c) were not found in the Na records. Ac- 
cording to Sodemann and Stohl (2009), precipitation in this 
high-altitude region of the Eastern Antarctic originates from 
sources located much further north than for the coastal re- 
gions. The lack of a correlation between the annual rBC data 
at the remaining sites and the non-linear trends at 08-4, 08-5 
and 07-1 suggests that the low frequency rBC variability may 
be linked to emission variability, or site specific atmospheric 
transport. 

To test this hypothesis, the spectral coherence for Na and 
rBC was investigated (Fig. 6). This analysis determined the 
coherence between periodic signals in the rBC and Na time 
series. A high coefficient for a given frequency suggests that 
the two periodic signals have coherent variability. Sites 08- 
4 and 08-5 exhibit coherence coefficients higher than 0.38 
(black line) for a large portion of the bandwidth. Notable 
exceptions are the ENSO band from ~4 to 7 yr and at lower 
decadal frequencies. For sites 07-2 and 07-7, coherence is 
much lower and often <0.38, confirming observations made 
on cross-correlations between Na and rBC (Fig. 5). Coher- 
ences between Na and rBC for sites 07-1 and 07-5, were sim- 
ilar to 08-4 and 08-5 with less coherence at low frequencies 
(<0.2 cycles per year). 

The spectral power of rBC time series is shown as a red 
line in Fig. 6. Peaks of high power designate frequencies 
explaining some variability of the signal. If these power 



Frequency (cpy) 


o 

u 


Fig. 6. Spectral power (red) of rBC NUS records and coherence 
coefficients (black) between rBC and Na investigated for the whole 
period (since 1800). Non-zero coherence is above 0.38. Red letters 
“NT” stand for “Non Transport”. They indicate periodic signals in 
the rBC records that are not coherent to Na (no black peak) and 
that are likely related to rBC emissions rather than regional to long 
range atmospheric transport. The red numbers below "NT” show 
the corresponding periodicity (in years). 


peaks do not correspond to a peak in coherence between Na 
and rBC (black line), they indicate an oscillation that is not 
linked to common atmospheric transport (NT). A periodicity 
of ~4.5 to 7 yr (~0. 17 to 0.2) is found common to sites 07-2, 
07-5, 07-7, 08-4 and 08-5. This period window suggests the 
influence of ENSO. However, even if an ENSO “signature” 
is present, none of the rBC records is statistically correlated 
to the ENSO index, which may be explained by two reasons. 
First, the records do not have the temporal resolution to ade- 
quately resolve the signal from noise. Second, ENSO has ,by 
nature, a dual effect on fire potential, by inducing drought on 
one side of the Pacific and floods on the other side, rendering 
a potential ENSO-fire signal very disparate (Rrawchuk and 
Moritz, 2011). 

However, a NT periodic oscillation of ~15 to 40 yr 
(0.05 ± 0.024 cycles per year) is found in records 07-1, 07-2, 
07-5, 07-7 and 08-4 (Fig. 6). This long-term periodicity in 
rBC, which is not related to Na, suggests a link to fire emis- 
sion variability or long-range upper atmospheric transport. 
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However, both may be linked to ENSO. Indeed, this periodic- 
ity is found to correspond to an ENSO reconstruction derived 
from the North American Drought Atlas (NADA) (Li et al., 
2011), which is anti-correlated to the large-scale rBC vari- 
ability observed in Fig. 3d (scale inverted). Here, higher rBC 
concentrations are associated with low variance periods (La 
Nina, colder) and lower concentrations during high variance 
periods (El Nino, warmer). This suggests that increased rBC 
emissions from fire drive higher rBC loading of the Antarc- 
tic atmosphere during decadal time periods dominated by La 
Nina. 

4 Conclusions 

Concentrations of rBC found in the NUS ice cores reveal 
both spatial and temporal variability during the 1800-2000 
time period. Spatial variability was primarily associated with 
changes in elevation and is likely linked to increased atmo- 
spheric loading in rBC and/or decreased accumulation with 
altitude. Relatively stable net snow accumulations rates at 
the NUS sites (Anschutz et al., 2011) suggest that decadal 
variability is related to changes in the rBC aerosol in the 
overlying air. On the other hand, the absence of strong cor- 
relations between the records may suggest that site-specific 
atmospheric transport and surface processes influence rBC 
concentrations at these sites. This is confirmed by high 
correlation and coherence coefficients between Na and rBC 
for some of the sites. This observation is different from 
the results obtained at WAIS and Law Dome by Bisiaux 
at al. (2012). Indeed, at those low elevation sites, it was 
shown that most of the recorded rBC variability was inde- 
pendent from atmospheric transport, which modulates the Na 
record. However, spectral analysis revealed the existence of 
non-transport oscillatory signals common to almost all the 
records. Common features in the records’ non-linear trends, 
showing relatively low concentrations from 1890 to 1910, 
high concentrations until 1930 and an increasing trend at the 
end of the 21st century, confirm the presence of a variability 
linked to rBC sources only. Nevertheless, while large-scale 
changes in rBC deposition at WAIS and Law Dome were 
found to correspond to a change in anthropogenic activities, 
measurements from the East Antarctic Plateau suggest a link 
with ENSO long-term emissions. In any case, global climate 
and aerosols models may enlighten the variability of rBC de- 
position to Antarctica, and the apportionment between the 
various continental sources. 
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